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ABSTRACT 

Context. Recent reaction rate evaluations include reaction rate uncertainties that have been determined in a statistically meaningful 
manner. Furthennore, reaction rate probability density distributions have been determined and published in the form of lognormal 
parameters with the specific goal of pursuing Monte Carlo nucleosynthesis studies. 

Aims. To test and assess different methods of randomly sampling over reaction rate probability densities and to determine the most 
accurate method for estimating elemental abundance uncertainties. 

Methods. Experimental Monte Carlo reaction rates are first computed for the -^Ne+a, ^''Ne(p,7)^'Na, -'Mg(p,7)^^Al, and '*'F(p,t!')'^0 
reactions, which are used to calculate reference nucleosynthesis yields for 16 nuclei affected by nucleosynthesis in massive stars and 
classical novae. Five different methods of randomly sampling over these reaction rate probability distributions are then developed, 
tested, and compared with the reference nucleosynthesis yields. 

Results. Given that the reaction rate probability density distributions can be described accurately with a lognormal distribution, 
Monte Carlo nucleosynthesis variations arising from the parametrised estimates for the reaction rate variations agree remarkably 
well with those obtained from the true rate samples. Most significantly, the most simple parametrisation agrees within just a few 
percent, meaning that Monte Carlo nucleosynthesis studies can be performed reliably using lognormal parametrisations of reaction 
rate probability density functions. 

Key words. Nuclear reactions, nucleosynthesis, abundances - Methods: numerical 



1. Introduction 

With the growing availability of high power computers, our abil- 
ity to explore the effects of nuclear physics uncertainties on nu- 
cleosynthesis has reached new levels. With post-processing cal- 
culations, for example, thousands of computations can be per- 
formed in relatively little time, thus opening new avenues for 
exploration. Furthermore, better treatment of reaction rate uncer- 
tainty propagation has recently lead to the publication of rates 
that include statistically meaningful uncertainties (Longland' 



|et al.j |2010[ [Tliadis et al., ,2010c_a.b) and include lognormal rep 



resentations of the rate probability density function^ Given that 
information, methods should be developed to utilise those reac- 
tion rate uncertainties and investigate their elfect on nucleosyn- 
thesis. 

Monte Carlo methods are one option for investigating the 
effects of reaction rate uncertainties on nucleosynthesis yields. 
These methods can either be used to investigate the effects of a 
single reaction, a small group of reactions, or they can be used 
to investigate the overall effect of all reaction rates in complex 
environments involving a great number of competing reactions. 
The general strategy in the latter case is to (i) compute a sam- 
ple rate for every uncertain reaction simultaneously, (ii) compute 
nucleosynthesis yields using these samples in a post-processing 
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' Note that there is no restriction that the rates obtained in llliadisl 
|et al.|p010c^ must be within the published Icr "high" and "low" rates 
but rather follow a continuous probability density distribution. 



model, and (iii) repeat many times. Correlations between abun 
dances and reaction rates can be found in a post-analysis proce 
dure to identify reaction rates of interest. These methods have 
be used previously in a number of studies (e.g., Hix et al.||2003 



Roberts , 2006 ; Par ikh et al.| |2008 ). In each of those studies, tem 
perature independent enhancement factors were applied to the 
rates of all nuclear processes to investigate their effect on nucle- 
osynthesis in novae and x-ray bursts. The enhancement factors 
were sampled from lognormal distributions that were chosen to 
represent the estimated uncertainty of the reaction rates at tem- 
peratures most impoitant in these environments. In some studies, 
these estimations were arrived at by considering very general 
characteristics of the reaction rates such as whether they are de- 
rived purely from theoretical predictions or were measured with 
radioactive beams (e.g., |Hix et aL] |2003[ ). Until now, however, 
these methods have not considered the rate uncertainties for each 
reaction separately. 

The method of varying reaction rates by a constant factor 
assumes that the tiTie, unknown, reaction rate differs from the 
recommended rate by a constant factor over all temperatures. 
Equivalently, this method makes the assumption that rate un- 
certainties are temperature-independent. However, in reality, the 
uncertainties of individual reaction rates are, indeed, often tem- 
perature dependent if they are based on experimental constraints, 
usually with larger uncertainties associated with the rate at lower 
temperatures. This can clearly be seen in the uncertainty bands 
presented by Iliadis et al. (2010b). Furthermore, by considering 
the way that reaction rates are computed, it is clear that the tem- 
perature dependence of a sample reaction rate should not nec- 



1 



R. Longland: Recommendations for Monte Carlo nucleosynthesis sampling (RN) 



essarily follow the same temperature dependence as the recom- 
mended rate. For example, the rate could be higher than the rec- 
ommended rate at some temperatures, while being lower at oth- 
ers. How important is this effect? The et al. (2000) discussed the 
"race against time" effect that different temperature dependen- 
cies for the ^^Ne+a reaction rates might have on the s-process in 
massive stars. Presumably, the complex interaction and compe- 
tition of reaction rates in a dynamically changing environment 
means that simple enhancement of reaction rates does not fully 
explore the nucleosynthesis variations arising from reaction rate 
uncertainties. More advanced methods should, therefore, be in- 
vestigated to fully account for this effect in nucleosynthesis un- 
certainty studies. The rate probability density distributions pub- 
lished in |Iliadis et al. (2010ci make this investigation possible, 
while eliminating the need for simplifying estimations of rate 
uncertainties based on general reaction properties. 

The obvious choice for correctly performing Monte Carlo 
nucleosynthesis studies is to directly use Monte Carlo reac- 
tion rate samples obtained separately from the RatesMC code 
(Longland et al. 2010| l. In this case, individual samples of the 
nuclear physics input are used, thus accounting for all possible 
behaviour of the reaction rates. However, this requires a consid- 
erable amount of effort using tools that most computational as- 
trophysicists do not have access to. The purpose of this research 
note, therefore, is to investigate different rate sampling schemes 
in comparison to the ideal case discussed above, and present the 
best practice with which to utilise the reaction rate uncertain- 
ties that are published in Iliadis et aL] ( |2010 c ). In Sec. |2] a num- 
ber of possible sampling schemes are discussed. Those schemes 
are tested with an example set of reaction rates and discussed in 
Sec.|4] Conclusions are given in Sec.|5] 



2. Reaction rate sampling metliods 



The reactions evaluated by Iliadis et al. ( 2010c[ ) involved target 
nuclei with masses in the A= 14 - 40 range, whose cross sec- 
tions are usually dominated by resonant capture. For reactions 
that involve a large number of resonances, it is clear that a great 
number of possibilities exist for how a sample reaction rate could 
vary with respect to its recommended value. In light of this, the 
most consistent way to sample the reaction rates is to sample the 
nuclear physics input separately for each nucleosynthesis calcu- 
lation. Since most computational investigations do not have ac- 
cess to the data necessary to perform this kind of sampling, this 
method is viewed as the ideal, yet unrealistically complicated 
standard. For the purpose of this study, the RatesMC code used 



by Longland et al. (2010 1 is used to compute this ideal case for 
comparison alongside the estimations considered in this work. 
Throughout the following discussion, rate variations obtained 
from samples of the nuclear physics uncertainties are referred 
to as the "optimum" rate samples. 

In order to construct rate variation schemes, we must first 
evaluate the expected behaviour of the reaction rates as a func- 
tion of temperature. A reaction rate per particle pair is computed 
by 



(crv) 
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where yU is the reduced mass of the reacting particles, ^ — 
M()Mi /(Mq+Mi); Mi denotes the masses of the particles; k is the 
Boltzmann constant; E is the centre-of-mass energy of the react- 
ing particles; and o-(E) is the reaction cross section at energy. 



E. This cross section frequently consists of a collection of reso- 
nances, so the reaction rate varies according to resonance prop- 
erties, such as resonance strength or energy. The convolution of 
these resonance cross sections with the Boltzmann distribution 
in Eq. ([T]) will always result in a smoothly varying reaction rate 
with temperature. 

To construct a rate variation scheme, first consider the avail- 
able information from Iliadis et al.| ( |2010cj ); information that is 
also included directly in the starlib library (see, for example, 
[Iliadis et"aL 2011 1. The tabulated rates include two parameters, 
fi and cr, to describe the rate probabiUty density distribution at 
each temperature (see Longland et al. 2010[ ): 
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where Xmed, JCWgh, and xiov, are the recommended (median), high, 
and low reaction rates, respectively. The variables // and cr here 
refer to the lognormal shape parameters, and are not to be con- 
fused with the more commonly used Gaussian mean and stan- 
dard deviation. A reaction rate sample, x{T), at a specific tem- 
perature, T, can be represented by 



x{T) = e^'^' > ■ e- 



.KT) . „p{T)a-{T) 



(3) 



where p{T) is a random variable that is normally distributed (i.e., 
distributed according to a Gaussian distribution with an expec- 
tation value of and standard deviation of 1). The second com- 
ponent of Eq. ([3]) is dubbed as the "uncertainty factor", and is 
a temperature dependent quantity. Given this information, the 
challenge becomes one of choosing values for p{T), recalling 
that these values must vary smoothly with temperature, and must 
be distributed normally to enforce the lognormality of reaction 
rate probability densities. 

The simplest parametrisation for p{T) is determined by as- 
suming that it is independent of temperature, i.e., pij) - a 
where a is sampled from a normal distribution. This is dubbed 
the "flat" parametrisation in the following discussion, although 
it is important to emphasise that for this parametrisation, the 
uncertainty factor in Eq. (|3]l is still temperature dependent, and 
given by e""'^^''. This was not the case in the Monte Carlo varia- 
tions of previous studies (e.g., |Parikh et al. 2008)). Those studies 
applied an uncertainty factor to the rate regardless of the un- 
certainty temperature dependence. Their method is equivalent to 
usinga constant value, /, for the uncertainty factor In this case, 
Eq. (Bb would become x{T) - f e'^'-^K As mentioned before, this 
is not representative of the rate uncertainty at all temperatures for 
experimentally constrained reaction rates. It is, however, a good 
approximation for reaction rates derived purely from theory. 

Consider, now, the case in which the reaction rate uncer- 
tainty can be characterised in two regions. A single rate sample 
could, for example, be high with respect to the recommended 
rate at lower temperatures, while being low at higher tempera- 
tures. This situation can be reproduced by parameterising p(T) 
in Eq. ^ with a hyperbolic tangent function. To further gener- 
alise this parametrisation, an offset parameter can be utilised: 



p{T) - lo + a tanh 

V2 



2.94 



logio(rjr) 



(4) 



where the sampled variables a, o, S , and are shown in the in- 
set of Fig. [T] and are defined as the maximum deviation, offset, 
spread parameter, and cross-over temperature, respectively. The 
spread parameter is defined as the logarithm of the temperature 
range needed for p to vary from p - o + Q.9a io p - o - 0.9a. 
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Fig. 1. (Colour online) Example of a single rate sample of the 
^^Ne(Qf,y)^^Mg reaction using the parametrised rate factor vari- 
ations. Plotted is the reaction rate ratio (i.e., compared to the 
recommended rate). Therefore, the line at unity represents the 
recommended rate, while a point at 10 means that the rate is 
10 times that of the recommended rate. The red density distri- 
bution represents the coverage probability of the reaction rate, 
with thick and thin black lines denoting the 68% and 95% un- 
certainties, respectively. The dashed blue line represents a single 
sample of the reaction rate obtained when using Eq. (|3]l with 
p(T) determined from the parametrisation in Eq. Q and shown 
in the inset. 



In this parametrisation scheme, a and o must be sampled from 
normal probability density distributions, while S and Tc can take 
on a range of values sampled from flat probability density func- 
tions. Here, the ranges are chosen to be 0.1 < S < 1 GK and 
10"^ < < 10 GK. This parametrisation is referred to as 
the "full" parametrisation hereafter. An example of applying this 
parametrisation to the ^'Ne(Q',n)^^Mg reaction rate is shown in 
Fig.[T] where a = -1.5, o = -0.6, = 0.1 GK, and 5 = 0.4 GK. 

To further investigate the hyperbolic tangent parametrisation 
of the reaction rate samples, two more cases are considered. 
In the first case (the "no offset" case), the offset parameter is 
constrained to zero. This requires removing the 1 / V2 factor in 
Eq. Q to ensure normality of the distribution. The parameter 
space for the reaction rate will be constrained for this case, in 
particular, the possibility that the rate has a systematic shift at 
all temperatures will not be accounted for. The nucleosynthesis 
yield uncertainties determined from this method are therefore 
expected to be reduced using this parametrisation. For the sec- 
ond case (called "fixed S and Tc"), the parameters, S and Tc 
are fixed at values of 5 = 0.2 and = 0.2 GK. These val- 
ues are chosen to maximise the reaction rate shape variations 
around 0.2 GK. This final case should help answer the question 
of whether the temperature dependence of the ^^Ne+a reactions 
has a large effect on s-process yields as predicted by The et al. 
( |2000l l. 



Finally, although we have argued that p(T) must vary 
smoothly with temperature, it can be illustrative to investigate 
the effect on nucleosynthesis when p(T) this constraint is re- 
laxed. On a grid of temperatures, therefore, p(T) is sampled 
independently so that the resulting rate can fluctuate randomly 
from one temperature to the next. These sample rates are ob- 
tained once per Monte Carlo nucleosynthesis computation. This 
will be referred to as the "random" case. Another technique is 
to sample all reaction rates at every evaluation step within the 
nucleosynthesis code. In this case, the rate will change rapidly 
around the recommended rate over small temperature ranges, 
thus yielding unrealistically small nucleosynthesis uncertainties. 
This latter case is not investigated here. 

3. Test Cases 

The effects of the 5 rate parametrisation schemes, in addition to 
the optimum rate samples calculated directly from the nuclear 
physics input, are investigated with respect to a test suite of re- 
actions and scenarios. To achieve this, post-processing calcula- 
tions are used to calculate Monte Carlo nucleosynthesis yields 
with some representative profiles and networks. To obtain a good 
representation of the uncertainties, 3000 Monte Carlo samples 
are used for each method, which are found to reproduce abun- 
dances to within 3% and uncertainties to within about 10% be- 
tween runs. 

The first of the test cases are the ^^NeH-a reactions (i.e., 
^^Ne(c!',y)^^Mg and ^^Ne(Q',n)^^Mg), which are important for s- 
process neutron production in massive stars and asymptotic gi- 
ant branch stars. For an overview of the reactions and their nu- 
clear physics uncertainties, the reader is referred to, for example, 
IJaeger et al.l ( [200T) ; [KoehIerl ( [2002) i; [Longland et al.| 



(20121. A 



temperature-density profile corresponding to core helium burn- 
ing in a 25 Mq massive star is used to study nucleosynthesis for 
these reactions, and has been discussed elsewhere by The et al. 



(2000 1 and Iliadis (2007j|. While helium burning occurs through- 
out the profile, the^^Ne+o- reactions are only active towards the 
end when temperature and density become high enough. A nu- 
cleosynthesis network containing 583 nuclei is used, which con- 
tains nuclei around the nuclear valley of stability from A = 1 to 
A - 100. This network is sufficient to study the weak compo- 
nent of the s-process. To characterise the nucleosynthesis yield 
uncertainties, final abundances of 4 nuclei are considered: ^^Mg, 
^^Cu, ^"Zn, and **^Rb. These are chosen to represent nuclei that 
are affected by the ^^Ne-Ha reactions. 

The other three reactions considered here are important in 
nova nucleosynthesis. We consider a single zone temperature- 
density profile from the hottest hydrogen-burning zone that cor- 
responds to a 1.25Mq ONe white dwarf accreting material at 
a rate of 2 x 10""'Mo.y~'. This model is similar to the "PI" 



model considered by Iliadis et al. (2002 1 with initial abundances 
equal to those presented in their Tab. 2. A nucleosynthesis net- 
work containing 146 nuclei from A = 1 to A = 50, and linked 
by 1283 nuclear processes is used to integrate the system's evo- 
lution during a single nova outburst. Two species of ^^Al are 
used to represent the abundances of ^^Al in its ground state, 
-^Al*, and its isomeric state, ^''Al'". It was shown previously 
by Iliadis et al. (2002) that some important reactions for this 
model are "'F(p,Q')'"0, 20Ne(p,r)2iNa, and 25Mg(p,7)2''Al. The 
first of these has large uncertainties arising from ambiguities in 
the nuclear physics data (see, for example Beer et aL]|2011[ ), and 
strongly affects the final yields of oxygen and fluorine isotopes. 
The second reaction considered is important because at the start 
of the nova explosion, ^"Ne comprises 25% of the material's ini- 
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tial abundance. The reaction proceeds relatively slowly, thus it 
strongly affects abundances throughout the neon-silicon region. 
The final reaction considered is important because it directly af- 
fects the abundances of magnesium and -^Al-''. These are impor- 
tant species in nova nucleosynthesis because they can be mea- 
sured precisely in pre-solar grains. The uncertainties for these 
reactions are adopted from |Iliadis et al.|p010a[ ). 

4. Results and discussion 

To illustrate the final distribution of abundances that arise from 
the Monte Carlo nucleosynthesis calculations, the final abun- 
dances for ^^Mg obtained from the optimum rate samples, the 
full and flat parametiisation schemes, and random sample meth- 
ods are plotted using cumulative frequency plots in Fig. |2] A 
number of features are immediately obvious. First, the median 
abundances (i.e., the point at which the cumulative distribution 
crosses the 50% point) agree within statistical fluctuation for all 
methods, as is expected. This is because, provided proper sam- 
pling has been carried out, the median abundances will corre- 
spond to the median rate regardless of the parametrisation used. 
Secondly, the random parametrisation produces a significantly 
wider distribution than the results from optimum rate samples. 
The unphysical, rapidly changing reaction rates possible when 
using this method could be responsible for this, thus producing 
hard-to-predict nucleosynthesis results. Indeed, we find that for 
other reactions such as '^F(p,Q;)'^0, the opposite effect occurs, 
i.e., the uncertainties resulting from this method are smaller than 
those obtained from optimum rate samples. We conclude, there- 
fore, that this random sampling of rates is not only unphysical, 
but also unreliable for calculating accurate nucleosynthesis yield 
uncertainties. Thirdly, full parametrisation of the reaction rate 
samples using the hyperbolic tangent function in Eq. Q agrees 
remarkably well with the results from using the optimum rate 
samples, albeit with slightly underestimated uncertainties. 

Perhaps even more remarkable from Fig.|2]is that the uncer- 
tainties derived by using flat distributions for p(T) also agree ex- 
cellently with the uncertainties derived from using optimum rate 
samples, and are barely distinguishable from the full parametri- 
sation method. This finding is further enforced in Fig.|3] where 
the uncertainty ranges obtained for key isotopes arising from 
varying the four reactions described in Sec. [3] are displayed as 
uncertainty ratios. To obtain these ratios, uncertainties are first 
determined from the 16th and 84th percentiles of their respec- 
tive cumulative distributions, and thus represent Icr nucleosyn- 
thesis yield uncertainties. These ranges are then normalised to 
the median abundance obtained from the optimum reaction rate 
samples. Only the ratios for optimum samples, full parametri- 
sation, random rate samples, and flat parametrisation are dis- 
played. The "fixed S and Tc" and "no offset" cases are not dis- 
played for reasons of clarity. For these two cases, the uncertainty 
ratios bands are only slightly smaller than those obtained from 
the full parametrisation for all cases considered here. This fig- 
ure confirms our findings that the full and flat parametrisation 
schemes reproduce the uncertainties obtained from optimum rate 
samples. The determined uncertainties are all within a few per- 
cent of those obtained from optimum rate samples, and certainly 
differ by smaller amounts than the uncertainties, themselves. The 
random parametrisation scheme consistently performs poorly, 
although for some nuclei (i.e., those not directly connected to 
the reaction), the effect is diluted somewhat. 

The only exception to these findings is the '*'F(p,Q')'^0 re- 
action, for which the uncertainties are too large by roughly a 
factor of two for "^O and '^F when utilising the parametrisa- 
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Fig. 2. (Colour online) Final ^^Mg abundance cumulative den- 
sity plots the three methods: (i) Input of reaction rate samples 
obtained directly from RatesMC ("Optimum Samples"), (ii) full 
parametrisation of p from Eq. Q ("Full"), (iii) randomly sam- 
pled reaction rates at each temperature ("Random"), and (iv) a 
constant value for as a function of temperature ("Flat"). 



tion schemes. The random parametrisation scheme, contrary to 
the other reactions, yields smaller nucleosynthesis uncertainties 
in this case. The '**F(p,ci')'^0 reaction presents a challenge be- 
cause the cmTent uncertainty in the rate is dominated by un- 
known factors in the interference effects between a number of 
2010a| l. The reaction rate probability 



resonances (Iliadis et al. 



density distribution is, therefore, bi-modal in shape, and hence 
poorly represented by a lognormal distribution. This fact is re- 
flected in the poor Anderson-Darling statistics at most relevant 
temperatures for this reaction. The Anderson-Darling statistic 
represents a measure of how well the true reaction rate distribu- 
tion is reproduced by the lognormal approximation, and is dis- 
cussed in detail in Longland et al. ( 2010| l. In this case, therefore, 
the cr parameter describing the rate uncertainty is a bad approx- 
imation, so using it to compute reaction rate uncertainties is ex- 
pected to fail in reproducing the optimum sample uncertainties. 
However, for this case, the poor lognormal shape of the reac- 
tion rate distribution already prompts further experimental in- 
vestigation. With the addition of new, more precise nuclear data, 
the rate will begin to approach a lognormal distribution more 
closely, and hence the parametrised Monte Carlo approach for 
this reaction will perform better. 

Finally, the smallest abundance uncertainties for all nu- 
clei considered come from using the fixed value for S and 
parametrisation. This arises from the values adopted for the fixed 
parameters, which were chosen to maximise the reaction rate 
shape effect discussed before. Although the shape of the rates 
will be changing as a function of temperature, the magnitude of 
the rates will, on average, be closer to the recommended rate 
around T= 0.2 GK. We can conclude from this that, at least for 
the present cases, the shape of the reaction rate has little effect 
on nucleosynthesis yields. 
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Fig. 3. (colour online) The uncertainties computed from 3 of the parametrisation schemes investigated here compared to the individ- 
ually sampled reaction rates. The final abundances for four nuclei are considered for each reaction under investigation here. Each 
group of vertical bars for each nucleus represents, from left to right, the abundances calculated using: (i) The "optimum" reaction 
rate samples; (ii) the full parametrisation; (iii) the unphysical, random, case; and (iv) constant values of p(T) - a. 



5. Conclusions 

Reaction rates provided by the recent reaction rate evaluation 
of Iliadis et al. ( 2010c| l included temperature dependent lognor- 
mal parameters that were found to accurately reflect the reac- 
tion rate probability densities. However, to date, no investiga- 
tion has been presented that determines the most effective way 
to use these parameters. Past research (The et al. 2000| l showed 
that nucleosynthesis yields can be influenced not only by reac- 
tion rate magnitudes, but also by their temperature dependence. 
An investigation into how best to randomly sample the reaction 
rates provided by Iliadis et al. ( 2010c| l is therefore performed 
here. The investigation was undertaken for a range of reactions 
occurring in core helium burning in massive stars (^^Ne-no') and 
in novae (20Ne(p,y)2iNa, 25Mg(p,r)2<'Al, and ^^F(p,a)^^0). 

By applying a temperature-dependent rate multiplication 
factor, e''(^)°'(^), to the reaction rates, the isotopic yields follow- 
ing post-processing of massive star helium burning and nova ex- 
plosion models were investigated for a range of p{T) parametri- 
sation schemes. These choices were compared with the yield un- 
certainties obtained when reaction rate samples are obtained di- 
rectly from the nuclear physics input and calculated using the 
RatesMC reaction rate code from Longland et al. ( 2010| l. 

The most significant finding was that not only can the op- 
timum reaction rate sample uncertainties be reproduced well 
by using a hyperbolic tangent function for p{T), but they are 
approximated remarkably well by a flat function: p{T) = a. 
The temperature dependence of the reaction rates investigated 
were, therefore, found to be relatively unimportant to nucleosyn- 
thesis yields. This finding is particularly significant because it 
means that when performing nucleosynthesis sensitivity studies, 
all reaction rates can be sampled simultaneously. The sensitivity 
of nucleosynthesis yields to particular reactions can be charac- 
terised by plotting abundance yields against p (in a similar fash- 



ion to Fig. 8 in Parikh et al. (2008 1) and computing measures of 
correlation such as the Pearson or Spearman correlation coeffi- 
cients. 

For investigating the efl'ect of single reaction rate uncertain- 
ties, or small groups of reactions where coiTelation coefficients 
are not needed, adopting either the full, or flat parametrisation 
schemes is recommended. Using the full parametrisation scheme 
is useful for investigating the effect that temperature dependence 
of a reaction rate has on nucleosynthesis. For rates with large 
uncertainties, such as those involving unstable target nuclei, this 
could still have an effect, and should always be considered when 
evaluating reaction rate uncertainties. 
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